Democratic particle motion for meta-basin transitions in simple glass-formers 
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We use molecular dynamics computer simulations to investigate the local motion of the particles in 
a supercooled simple liquid. Using the concept of the distance matrix we find that the a— relaxation 
corresponds to a small number of crossings from one meta-basin to a neighboring one. Each crossing 
is very rapid and involves the collective motion of O(40) particles that form a relatively compact 
cluster, whereas string-like motions seem not to be relevant for these transitions. These compact 
clusters are thus candidates for the cooperatively rearranging regions proposed long times ago by 
Adam and Gibbs. 
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In recent years significant progress has been made 
in our understanding of the relaxation dynamics of 
glass-forming liquids at intermediate and low tempera- 
tures. Sophisticated experiments and computer simu- 
lations have identified many of the salient features of 
this dynamics, and theoretical approaches have helped 
to rationalize them, at least to some extend 0, 0j IS Gi • 
Despite this progress, many of the most elementary 
questions have not been answered so far and among 
them is the nature of the motion of the particles in 
the a— relaxation regime at low temperatures. Experi- 
ments and simulations have demonstrated that this dy- 
namics is quite heterogeneous and therefore can be used 
to explain the observed stretc hing of the time correla- 
tion functions H, H |1 11 H ES El This hetero gencous 
dynamics has been shown to be related to cooperative 
motion in which a small number of particles (a few per- 
cent) undergo a collective relaxation dynamics in that 
they move, often in a string-like fashion, by a distance 
that is comparable to the one between neighboring par- 
ticles 0, 0, 0]. Since a qualitatively similar hetero- 
geneous dynamics has also been found in simple lattice 
models that show a glassy dynamics and for which it is 
well known that the a— relaxation is intimately connected 
to the dynamical heterogeneities (DH), it has been pro- 
posed that these DH are crucial for the relaxation dy- 
namics of all glass- forming systems 0]. However, since 
in these lattice models all elastic or quasi-elastic effects 
are completely neglected, it is not at all clear whether or 
not the DH are indeed the only relevant mechanism for 
the relaxation. 

Another approach to describe the relaxation dynamics 
on the time scale of the a— relaxation is by means of the 
so-called potential energy landscape (PEL) 0, 
(or more precisely the free energy landscape) and hence 
to describe the dynamics of the system by considering 
its trajectory in configuration space. This PEL is rugged 
due to the presence of barriers in the free energy and 



hence at low temperatures the resulting dynamics will 
be slow. Using the concept of the inherent structures, 
evidence has been given that (roughly speaking) the mo- 
tion of the system in the PEL can be decomposed into 
two types of movements 0,0: In the first type the sys- 
tem explores some minima which are locally connected 
to each other and are not separated by a significant bar- 
rier. Therefore such a collection of minima is called 
"meta-basin" (MB) and its exploration corresponds to 
the (3— relaxation. As a second type of motion, the sys- 
tem overcomes the barriers that surround a MB and en- 
ters a new MB, a motion that has been believed to corre- 
spond to the a— relaxation Q or to the elementary events 
of the a— relaxation [Tsj] . 

Although this picture for the motion of the system 
within the PEL is certainly appealing from a qualitative 
point of view, it does not give any insight on the na- 
ture of the dynamics of the particles on the microscopic 
level. Furthermore it remains unclear whether the DH 
mentioned above have anything to do with the barrier- 
crossing of the system moving in the PEL. The goal of 
the present work is therefore to clarify this issue and thus 
to make an advancement in our understanding of the re- 
laxation dynamics of supercooled liquids. To this aim we 
have done molecular dynamics computer simulations of 
a simple glass-former in order to identify the presence of 
the MB and to investigate the nature of the motion of the 
particles during the transition from one MB to another. 
The so obtained results can then be compared with the 
heterogeneous motion of the particles in order to see to 
what extend the two motions are related to each other. 

The system considered is a binary mixture of Lennard- 
Jones (LJ) particles. In previous investigations it has 
been shown that this system shows many features of 
glass-forming liquids and can thus serve as a simple 
model for such liquids |19| . The interaction between two 
atoms of type A (80%) and B (20%) is given by V a p{r) = 
4e Q /3{K/3/r) 12 - {<J al3 /rf}, where a,/3 e {A, B}. The 
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FIG. 1: Distance matrix A 2 (t', t") of the system for T = 0.50. 
The gray level correspond to values of A 2 (t' , t") that are given 
to the right of the figure. 

LJ parameters used are caa — 1-0, oaa = 1-0, e^B = 1-5, 
cab = 0.8, €bb — 0.5, and (Jbb — 0.88. These interac- 
tions have been truncated and shifted at r cuto ff = 2.5<j a 0. 
In the following we will use (Jaa and eaa as units of 
length and energy, respectively, and measure time in 
units of (ma\ A /^ e AA) 1 ^ 2 ■ The equations of motion were 
solved for the NVE ensemble at a particle density of 1.2, 
using the velocity form of the Verlet algorithm with a 
time step of 0.02. All the presented results correspond to 
the situation in equilibrium. 

In order to identify the MBs we define the following 
"distance matrix" (DM) [2(|: 

1 N 

A 2 (t',0 = ^El r '(*')-r i (0| 2 , (1) 

i=i 

where r<(i) is the position of particle i at time t. Thus 
A 2 (t',t") gives the system averaged squared displace- 
ment of a particle in the time interval that starts at t' and 
ends at t" . Note that the time average of A 2 (t',t' + 9) 
over t' gives the r— average of G s {r,6), the self-part of 
the van Hove correlation function for time displacement 
6. The same is true if one averages over a very large 
system. Since we are interested in individual MB-MB 
transitions, one has to avoid that the presence of several 
independent local rearrangements, that will occur in a 
large system almost simultaneously, obscures the analysis 
of the individual event. Therefore we have considered a 
rather small system of 150 particles, the smallest possible 
system that does not affect the interactions, i.e. the box 
size was two times r cuto ff for the A — A interaction [2l| . 

Figure^shows a typical graph of the DM at T = 0.50 
as a function of the two time arguments t' and t" , with 
darker areas corresponding to configurations that have a 
smaller distance. At this temperature the dynamics is 
already slow with an a— relaxation time r of the order 
of 4 • 10 3 ( T can : e.g. be defined by requiring that 
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FIG. 2: Solid line (right scale): Averaged squared displace- 
ment S 2 (t, 6) for the trajectory given in 

Fig- El The value of 
is 160. Vertical bars (left scale): The function m(t,0) which 
gives the fraction of mobile particles, i.e. particles that moved 
more than the threshold value r t h = 0.3 in the time interval 
[t,t + 0], using 6 = 40. 

the self intermediate scattering function F s (q, t) has de- 
cayed to 10% of its initial value.) From this figure we 
see immediately that the dynamics of the system is quite 
heterogeneous in time in that it stays for a significant 
time relatively close to one region in configuration space, 
dark square-like regions, before it finds a pathway to a 
new region. Thus this is clear evidence that the system 
explores the present MB before it moves on to a neigh- 
boring one. At this temperature the typical sojourn time 
within one MB is around 300-800 time units, which is 
around 10% of r. Thus this sojourn time corresponds 
to the time scale of t*, the time that previous investiga- 
tions have shown to be relevant for the dynamical hetero- 
geneities in the system and which is defined as the time 
at which one observes the maximum in the non-gaussian 
parameter a.2(t) la. |23l. and which at this temperature 
is 400 time units [ljj. Hence we can conclude that the 
MB-MB transitions arc relevant for the DH whereas it 
takes about 5-10 such transitions in order to make an 
a— relaxation. We also point out that from Fig. it be- 
comes evident that the time for a MB-MB transition is 
quite short, on the order of 100 time units, which thus 
corresponds to about 20% of t*. Thus the transition is 
significantly faster than the a— relaxation times r. 

In Fig. [2] we show for the same run and time inter- 
val 5 2 (t, 9), the (particle) averaged squared displacement 
(ASD) of the particles within a time interval 9 (solid 
curve). This function is defined as 

S 2 (t,9) = A 2 {t-6/2,t + 6/2) (2) 
1 N 

- ^Ei r *(*-0/ 2 )- r *(*+0/ 2 )i 2 ■ ( 3 ) 

i=l 

Thus 5 2 (t,9) is A 2 (t',t") measured along the diagonal 
t" = t' + 9 and hence the average of this quantity over 



3 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

r 



FIG. 3: The distribution function G s (r,t,t + 9) for different 
values of t (curves with symbols). The value of 9 is 40. The 
bold curve is G s (r,9), the self part of the van Hove function, 
a) Values of t in which the system is inside a MB. b) Values 
of t for which the system is about to leave a MB. 



different start times t gives the usual mean-squared dis- 
placement for time lag 9. A comparison of this ASD with 
Fig- shows that S 2 (t,9) is showing pronounced peaks 
exactly then when the system leaves a MB. Thus we see 
that changing the MB is indeed associated with a rapid 
motion as measured in the ASD. 

In order to understand the nature of the motion of the 
particles when the system leaves a MB, we have calcu- 
lated G 3 (r, t,t + 9), the distribution of the displacement r 
of a particle for a given time difference 9 = 40. (Note that 
the average of G s (r, t,t + 9) over t gives G s (r, 9).) This 
distribution is shown in Fig. for starting values t that 
correspond to times at which the ASD shows a plateau, 
i.e. when the system explores a MB. Also included in the 
graph is the self part of the van Hove function, G s (r, 9), 
and we see that the distributions G s {r, t,t + 9) are nar- 
rower and more peaked than G s (r, 9) thus showing that 
in a MB the system moves more slowly than on aver- 
age. In Fig. |2t> we show the same distributions but now 
for times in which the system is about to leave a MB 
(compare the values of t with Fig. For these values 
of t the distributions are displaced to the right with re- 
spect to G s (r, 9) , showing that in this time regime the 
motion of the system is faster than on average. Most 
noteworthy is the observation that this shift is relatively 



uniform, i.e. a substantial part of the particles moves 
quicker than on average. Thus we can conclude that the 
rapid increase of the ASD is not due to the presence of a 
few fast moving particles, but instead to a "democratic" 
movement of many particles, in contrast to the results 
for cooperative motion on the time scale of t* which has 
been documented in earlier work HE 03. Thus this 
movement is very different in nature from the "string- 
like" motion found in the context of the dynamical het- 
erogeneities [l2ll23| . 

To demonstrate that the number of particles that par- 
ticipate at this democratic motion is indeed substantial 
and strongly correlated with a strong increase in the 
ASD, we have defined as "mobile" all those particles 
that in the time interval 9 = 40 have moved more than 
r t h = 0.3, and denote the fraction of such particles by 
m(t,6) [24J- I n Fig. [21 we have included the fraction of 
mobile particles as a function of time (vertical bars) and a 
comparison of this data with the ASD in the same graph 
shows that the fraction of mobile particles is indeed large 
whenever the ASD increases rapidly. This fraction is on 
the order of 30% of the particles and thus significantly 
larger than one would expect from G s (r, 9) if one inte- 
grates this distribution from r t h to infinity and which 
gives 0.09. 

In order to give an idea on the nature of the motion 
of the mobile particles during a MB-MB transition, we 
show in Fig. 0] a typical configuration of mobile particles 
before such a transition event and attach to each particle 
an arrow which points to the location of the particle after 
the event, i.e. a time 9 = 40 later. From this graph we 
recognize that the MB-MB transitions correspond to a 
movement in which the particles form a relatively com- 
pact cluster. Thus this is in contrast to the type of mo- 
tion found in the context of the DH in which the particles 
form string- like objects [T3, l25| . These compact regions 
have, at the temperature considered, around 30-60 parti- 
cles and can be considered as potential candidates for the 
cooperatively rearranging regions (CRR) proposed long 
time ago by Adam and Gibbs |26( and which are also at 
the heart of the approach of Goldstein for the relaxation 
dynamics of glass- forming systems [27^. 

Finally we mention that we have found that upon a 
decrease of temperature the sojourn time of the system 
within one MB increases rapidly. This is in agreement 
with previous results in which the concept of inherent 
structures was used ^3 , although here we have used 
a significantly simpler method to identify the MB and 
which can notably also be used in real experiments such 
as colloidal systems 0, E3- (We emphasize, however, 
that we have obtained qualitatively the same results by 
considering the inherent structures, although this ap- 
proach is computationally much more involved.) On the 
other hand an increase of the temperature makes that 
the structure of the MB is basically washed out and the 
ASD does no longer show the pronounced peaks. 

The present result indicate that the a— relaxation is 
not directly related to the presence of strings that are 
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FIG. 4: Configuration snapshot of mobile particles occurring 
in the MB-MB transition t = 680 -> t = 720. The spheres 
(light and dark for the A and B particles, respectively) give 
the location of the particles before the rearrangement and the 
arrows point to their position after the transition. The black 
spheres give the location of the corner of the simulation box. 

formed by a small number of particles but instead is due 
to a cooperative rearrangement of a substantial fraction 
of the particles. This cooperative rearrangement is re- 



sponsible for the transition between adjacent MBs and 
involves, at the temperature considered, on the order of 
40 particles. During such a MB-MB transition the major- 
ity of the particles outside this CRR does not contribute 
significantly to the relaxation. Therefore one needs on 
the order of 5-10 such transition events in order to com- 
plete an a— relaxation on the scale of the local neighbor- 
hood. Since these transitions take only a few percent of 
the a— relaxation time, one thus can envision, at least 
at low temperatures, the a— relaxation as a sequence of 
rapid, localized, cooperative relaxation events in which a 
substantial number of particles participate. This result 
is thus in surprisingly good agreement with the picture 
put forward long time ago by Adam and Gibbs and Gold- 
stein. 

Finally we mention that in view of the proposed con- 
nection between the geometrical properties (distribution 
of barrier heights, size of meta-basins,...) and the fragility 
of glass- forming systems 0, , it will be of great inter- 
est to do the analysis presented here for a strong glass 
former. This will allow to obtain a better understand- 
ing of the difference between strong and fragile glass- 
formers [2^| . 

GAA is research fellow of CONICET and RAM is re- 
search fellow of CIC. Financial support from Fundacion 
Antorchas, CONICET, UNS, SeCyT, DYGLAGEMEM, 
and ECOS is gratefully acknowledged. 



[1] M. D. Ediger, C. A. Angell, and S. R. Nagel, J. Phys. 

Chem. B 100, 13200 (1996). 
[2] W. Gotze, J. Phys.: Condens. Matter 10, Al (1999). 
[3] P. G. Debenedetti and F. H. Stihmger, Nature 410, 259 

(2001). 

[4] K. L. Ngai (Ed.): Proceedings of the Forth International 
Discussion Meeting on Relaxations in Complex Systems, 
Non-Cryst. Solids 307-310 (2002). 

[5] K. Schmidt-Rohr and H. W. Spiess, Phys. Rev. Lett. 66, 
3020 (1991). 

[6] S. Butler and P. Harrowell, J. Chem. Phys. 95, 4454 
(1991). 

[7] M. T. Cicerone, F. R. Blackburn, and M. D. Ediger, J. 

Chem. Phys. 102, 471 (1995). 
[8] W. Kob, C. Donati, S. J. Plimpton, P. H. Poole, and S. 

C. Glotzer, Phys. Rev. Lett. 79, 2827 (1997). 
[9] W. K. Kegel and A. van Blaaderen, Science 287, 290 

(2000). 

[10] E. R. Weeks, J. C. Crocker, A. C. Levitt, A. Schofield, 
and D. A. Weitz, Science 287, 627 (2000). 

[11] R. Richert, J. Phys.: Condens. Matter 14, R703 (2002). 

[12] C. Donati, J. F. Douglas, W. Kob, S. J. Plimpton, P. 
H. Poole, and S. C. Glotzer, Phys. Rev. Lett. 80, 2338 
(1998). 

[13] J. P. Garrahan and D. Chandler, Phys. Rev. Lett. 89, 

035704 (2002). 
[14] C. A. Angell, J. Non-Cryst. Solids 131, 13 (1991). 
[15] P. G. Debenedetti, T. M. Truskett, and C. P. Lewis, Adv. 

Chem. Eng. 28, 21 (2001); F. Sciortino, J. Stat. Mech. 

(in press) (2005). 
[16] C. Monthus and J.-P. Bouchaud, J. Phys. A: Math. Gen. 



29, 3847 (1996). 
[17] B. Doliwa and A. Heuer, Phys. Rev. E 67, 031506 (2003). 
[18] M. Vogel, B. Doliwa, A. Heuer, and S. C. Glotzer, J. 

Chem. Phys. 120, 4404 (2004). 
[19] W. Kob and H. C. Andersen, Phys. Rev. E 51, 4626 

(1995). 

[20] I. Ohmine, J. Phys. Chem. 99, 6765 (1995). 

[21] In order to check for the presence of finite size effects we 
have also repeated the presented analysis for a system 
with 1000 particles and considered only a subbox in this 
extended system. The results obtained were very simi- 
lar to the ones presented here and therefore finite size 
effects can be excluded. More details will be presented 
elsewhere [22^ . 

[22] G. A. Appignanesi, J. A. Rodriguez Fris, R. A. Montani, 
and W. Kob, to be published. 

[23] G. A. Appignanesi and R. A. Montani, J. Non-Cryst. 
Solids 337, 109 (2004). 

[24] We mention that if this threshold is changed by a rea- 
sonable amount, the number of mobile particles changes 
but the overall conclusions remain nevertheless valid. 

[25] G. A. Appignanesi, M. A. Frechero, L. M. Alarcon, J. A. 
R. Fris, and R. A. Montani, Physica A 339, 469 (2004). 

[26] G. Adam and J. H. Gibbs, J. Chem. Phys. 43, 139 (1965). 

[27] M. Goldstein, J. Chem. Phys. 51, 3728 (1969). 

[28] C. A. Angell, P. H. Poole, and J. Shao, Nuovo Cimento 
D 16, 993 (1994). 

[29] C. A. Angell, p. 1 in: K. L. Ngai and G. B. Wright (Eds.) 
Relaxation in Complex Systems (US Dept. Commerce, 
Springfield, 1985). 



